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i.  overview  and  observations 


Although  the  phenomenon  of  beam  plasma  discharge 
(BPD)  was  experimentally  observed  20  years  ago  l'2?  it  has 
recently  received  particular  attention  in  view  of  its 
importance  to  vehicle  neutralization  and  the  physical 
phenomena  exacted  in  ionospheric  electron  beam  injection 
experiments  from  rockets  and  the  space  shuttle  3-6.  while 
the  early  experiments1*2  injected  pulsed  keV  beams  into 
high  neutral  pressure  (p^l  mT)  high  magnetic  field  (.2-1=. 
KG)  devices,  more  recent  laboratory  experiments?”"11 
performed  at  the  large  tank  at  Johnson  Space  Center  (JSC), 
injected  steady  state  keV  beams  in  ionospheric  like 
pressures  (p^l-100  uT)  and  magnetic  fields  (B^.5-1.5  G). 
An  empirical  scaling  law  was  found  for  the  threshold 
current  Ic  required  for  ignition  of  BPD  as  shown  in  Fig.  1 
and  given  by  3/2 

Ic  S  ~x -  ftp'  (l) 

B  L 

where  is  the  energy  of  the  beam  and  L  the  length  of  the 
chamber.  The  value  of  A^.5-1  and  the  function  f (p)  has  a 
minimum  at  a  pressure  value  p0%l5yT  and  varied  roughly  as 
p^*5  to  p—1  above  and  below  p©  (Fig.  1).  This  experiment 
coupled  with  our  improved  understanding  of  the  nonlinear 
physics  of  beam  plasma  interactions  in  the  strongly 
turbulent  regime  12-15  allowed  us  to  make  major  progress 
in  the  physics  of  BPD.  The  interdisciplinary  nature  of 
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the  phenomenon  makes  it  an  appropriate  subject  for  a 
comment,  whose  purpose  is  to  delineate  the  plasma  physics 
of  the  B  PD  and  provide  a  simple  model  for  the 
incorporation  of  the  relevant  atomic  physics  and  transport 
processes. 

To  orient  the  reader  we  describe  first  the 
observations  of  BPD  at  the  JSC  chamber,  since  they  are 
rather  typical  of  similar  experiments.  An  energetic 
electron  beam  (E^  .7-1.5  keV,  I%l-100mA)  was  injected 
parallel  to  an  ambient  magnetic  field  (B^.8-1.5  G)  in  a 
chamber  of  axial  length  L  %20  m  and  approximately  equal 
diameter,  filled  with  neutral  gas  with  pressure  p%l-30uT. 
For  beam  currents  below  a  critical  value  lc»  which  is  a 
function  of  p,  L,  B  and  Eb  given  by  Equ.  (1)  the  beam,  as 
seen  optically,  follows  the  scalloped  trajectory  expected 
from  single  particle  dynamics,  accompanied  by  the  emission 
of  weak  whistler  and  cyclotron  waves.  For  beam  currents 
above  Ic  there  are  many  phenomena  indicative  of  strong 
collective  plasma  interactions,  such  as: 

(i)  The  beam  expands  radially  filling  its  radius  while  any 
sign  of  single  particle  dynamics  disappears. 

(ii)  Photometric  measurements  indicate  a  large  increase 
in  the  3914-A  volume  emission,  accompanied  by  jumps  in  the 
plasma  density  by  a  factor  between  5-10. 

(iii)  Electron  heating  to  1-2  eV,  accompanied  by  large 
fluxes  of  suprathermal  electrons  in  the  region  between  5- 
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(iv)  Large  localized  electrostatic  waves  near  the  plasma 
frequency  (u>e)  confined  inside  the  beam  regionr  while 
broadband  featureless  whistler  emissions  are  measured 
outside. 


II.  A  Simple  Physical  ££Q  Model: 


On  a  superficial  level  BPD  can  be  thought  as  anRF 
discharge  at  the  critical  density  (i.e.  “rf«“e)  but  with 
the  rf  fields  being  electrostatic  in  nature  and  driven  by 
a  beam  plasma  instability  (BPI).  In  probing  deeper  into 
the  BPD  observables  we  discover  many  important 
characteristics  different  from  RP  discharges.  Some  of  the 
key  plasmaphysics  issues  of  BPD  can  be  identified  by 
refering  to  the  diagram  of  Fig.  2.  In  the  absence  of  any 
collective  effects  the  system  will  follow  the  path  marked 
by©-  ©-  GX  Namely  the  beam  will  ionize  the  neutral  gas 
and  a  steady  state  will  be  established  by  balancing  the 
rate  of  ionization  with  the  rate  of  plasma  loss,  i.e. 


dn  „ 
dt 


n,  N<oV 

D 


n 

T 


(2) 


where  nb,  n,  N  are  the  beam,  plasma  and  neutral  density,  o 
the  ionization  cross  section  and  x  the  confinement  time. 
Although  we  will  neglect  recombination  here,  it  can  be 
included  by  taking  where  B  is  the  recombination 

coefficient.  The  expected  steady  state  density  from 
equ.(2)  is  then 
no 

r£  *  N<avb>  x  (3) 

If  r\,f(n ),  then  nQ  will  be  given  by  the  solution  of  equ. 
(3)  for  n.  In  order  to  have  BPD  the  path  marked  in  Pig.  2 
by©  -  ©-©should  be  allowed  on  a  time  scale  faster  than 
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t,  and  the  rate  of  ionization  due  to® should  be  larger 
than  due  to®  If  we  assume  that  the  path  is 
allowed  Equ.  (2)  should  be  modified  to  read 


n  m<°v  >-? 


(4) 


dt  "  “b"  “  b  “s"  's 

where  ns,  Vs  are  the  density  and  thermal  speed  of 
electrons  with  energy  above  the  gas  ionization  energy  E^. 
Assuming  that  the  confinement  time  is  not  affected,  the 


steady  state  density  after  BPI  will  be 

ll=  o 


n  <oV  > 

n,=  n_  (1+  ) 


v°v 


(5) 


n  <oV  > 

If  the  process  can  account  for  observations 

b  b 

(i),  (iii)  and  (iv)  above.  Namely  we  will  have  a  BPI  but 

n  <oV 

not  a  BPD.  If  — — rj^->l,  then  observation  (ii)  can  be 

n,<ov.  > 
b  b 

accounted  for  and  we  will  have  a  BPD.  The  above 
considerations  can  be  easily  extended  to  include  the 
possible  effects  of  BPI  to  the  plasma  confinement  (e.g. 
electron  heating)  by  following  the  path  marked  in 

Fig.  2.  Based  on  the  above  considerations  we  can  state  a 
generic  criterion  for  BPD  triggering.  "JBED  Mill  kfi 
triggered  ii  the  value  pi  the  px.Sr.fiED  density  n£  combined 
with  the  other  parameters  ( i.e.  p^j.  ££*  lu.  fix  fiLL  is  SUCh 
that  Ip  can  produce  pp  a  time  scale  fih,Qit.ei  than  lA  JuLULX 

pi  energetic  elections  such  that 


n  <oV  > 

-J _ £_  > 

nb<aVb> 


i: 
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It  is  also  easy  to  see  that  if  the  plasma  confinement 
during  BPI  is  not  affected  or  it  improves  we  can  reach 
stationary  state.  However  if  the  confinement  deteriorates 
the  final  state  will  have  relaxation  oscillations  with 
time  scale  the  new  confinement  time.  Therefore  the  plasma 
physics  input  to  the  problem  involves  the  computation  of 
conditions  for  BPI,  the  energy  transfer  efficiency  from 
the  beam  to  electrons  with  ionizing  energy,  and  the  plasma 
confinement  time. 
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III.  j!££  Scaling  BFD.L  The 


Guided  by  the  measured  plasma  physics  quantities  near 
BPD  threshold  at  the  JSC  experiment,  Papadopoulos4 
conjectured  that  the  BPD  threshold  is  associated  "  with 
the  value  of  the  plasma  density  n0  required  to  produce  an 
akaalii.fc.fi  instability  near  u»e%  The  reasons  for  the 
conjecture  were  the  following: 

(i)  The  value  of  the  electron  neutral  collisions  was  large 
enough  to  stabilize  kinetic  BPI.  The  hydrodynamic  BPI 
can,  however,  be  excited  due  to  its  larger  growth  as  well 
as  the  fact  that  it  is  a  negative  energy  wave  instability. 

(ii)  Convective  hydrodynamic  instability  has  group 


velocity  vg^Vb* 


sense 


that 


Since  the  system  size  L  is  short, 
^6-8,  and  the  instability  growth 


in  the 

Y<<  u)e 


convective  amplification  cannot  produce  substantial  energy 


deposition.  Only  an  absolute  instability  can  do  it. 


(iii)  The  reason  for  insisting  on  an  instability  near  u>e, 
is  not  only  the  experimental  evidence,  but  also  the 
requirement  for  fast  (collisionless)  energy  transfer  of 
the  energy  to  suprathermal  particles.  The  spiky 
turbulence  associated  with  we  waves  automatically  provides 
for  this.l*“15 


(iv)  An  assessment  of  the  overall  efficiency  of  energy 
deposition  by  the  beam  (i.e.  6-8%)  in  the  experiment  was 
consistent  with  stabilization  by  trapping. 


In  a  finite  size  system  an  instability  can  be  produced  by 
the  coupling  of  the  upper  or  lower  hybrid  plasma  wave  with 
the  negative  energy  (i.e.  slow)  wave  of  the  beam.  In  the 

weak  coupling  approximation  an  absolute  instability 
16 


requires 


v^v2  <  0 


,2>  V 


1 "  2 
1V2 


< !v,v,|) 


1/2 


Y 


=  L 


(6a) 

(6b) 

(6c) 


where  (vj_,  ^),  (V2 /  v2)  are  the  group  velocity  and 

damping  rate  (collisional  and  collisionless)  of  the  plasma 
and  slow  beam  wave.  The  first  criterion  comes  from  the 
requirement  that  the  unstable  wavepacket  encompasses  the 
origin  at  all  times.  The  second  is  the  need  to  overcome 
the  damping,  while  the  third  is  the  requirement  for 
positive  feedback.  Lc  corresponds  to  the  critical  length 
at  which  an  amplifier  becomes  an  oscillator  (i.e. 
oscillator  break  length).  Condition  (6a)  dictates  that 
the  instability  should  be  produced  at  the  upper  hybrid 
wave  (^ 0 )  which  is  a  backward  wave.  For  u0«ue  we  require 
-e>>ne.  The  growth  rate  (y)  and  vlf  v2  for  the  upper 

hybrid-beam  wave  intersection  are  qiven  by 

2  2 
R  cos  9 

D 

2 
o 

vr  vb  (7b) 

o2 

v2*  cos29sin29  (7c) 


1/2 


(7a) 
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where  R  is  a  geometric  fact or connected  with  the  finite 
radius  of  the  beam  (a)  and  the  plasma  (b),  wb  is  the  beam 
plasma  frequency  and 


2  _  i  ,  2  2. 

wo  2  (we  +  V 


[t<“ 


2  ,  n2,  202  2fl*l  1/2 

e  +  V^eV05  6J 


sin  0* 


kz  + 


(8a) 

(8b) 


The  value  of  ka  is  given  by  the  first  root  of  Jotk^b)*^, 
2  4 

i.e.  kx  = Criterion  (6a)  is  trivially  satisfied  for 
uQf  as  well  as  (6b);  the  condition  for  absolute 
instability  is  basically  (6c).  Prom  (6c)  and  (7)  we  find- 


b  ~  L 


(%\  1/2  P.eVb 

W  R1/2 


The  beam  density  is  given  by 


sin9 


(9) 


n, 

b 


a= 


I 

eVbTia2 

Vbsined 

?.e 


(10a) 


(10b) 


where  ^  is  the  accelerator  divergence  angle.  Noting  that 
in  the  experiment  a<b,  the  denisty  nQ  collisionally 
produced  by  the  beam  should  be  given  by  the  steady  state 
solution  for  a  line  source  of  the  diffusion  equationr  i.e. 


n  =  -122. 
o  2ueD 


(11) 


where  D  is  the  diffusion  coefficient.  In  comparing  Eqs. 

(11)  and  (3)/  we  find  that  they  are  equivalent  if  the 

a2 

confinement  time  is  diffusion  controlledf  i.e.  t  . 
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From  Eqs.  (9-11)  using  the  fact  that  o->-Eb,  we  find  the 


threshold  condition  as 


E3/2  D1'2 


sindsine 

where  K1  is  a  constant  including  the  term  ( - c - JL_), 

which  is  very  weakly  (i.e.  less  than  logarithmically) 
dependent  on  the  system  parameters.  It  was  shownlO  that 
in  the  low  pressure  regime  p<pQ  the  diffusion  is 
consistent  with  Bohm  diffusion  so  that 

d'db~5  p<Po  i13) 

Therefore  the  scaling  predicted  by  the  Papadopoulos 
conjecture  will  be 


C  p1/,2B'L^2L 


if  PQ  is  the  pressure  value  above  which  classical 
diffusion  dominates  (i.e.  DBs=D(p0))r  we  find  that 


D'Dc 


in  which  case  the  scaling  of  Ic  is 

E3/2  o 1/2 
Eb  ? 


Eqs.  (14)  and  (16)  are  consistent  with  the  experimental 
scaling.  In  fact  18  even  the  numerical  coeff icients,  which 
are  omitted  heref  accurately  reproduce  the  numerical 
values  of  lc.  Notice  that  the  relationship  Db»dcji  (Po) 
gives  po^B  which  seems  to  agree  with  Fig.  1. 


IV  .Concluding  Remarks 


We  proceed  next  to  discuss  the  applicability  of  the 
above  to  other  situations,  especially  the  ones  involving 
stronger  magnetic  fields.  We  first  notice  that  as 
discussed  before  p0%B;  we  therefore  expect  that  for 
situations  where  we>  increasing  the  magnetic  field  will 
increase  the  value  of  p0  at  which  the  transition  to  the 
high  density  regime  occurs.  Namely 

v  15  <§>  uT 

This  has  been  observed  by  Konradi  et  al19.  and 

Bernstein20,  where  for  B%3  8G  -r-  behavior  was  observed  up 

P' 

to  the  maximum  pressure  which  was  500  uT.  Notice  that  the 

observed  scaling  in  the  experiment  was  consistent  with 

— 5~ "5  '  as  given  by  our  conjecture19.  Similar  overall 
B*  P* 

scaling  was  found  in  Lyakhov  et  al21.  For  even  stronger 
magnetic  fields  when  U)e<i7e  at  threshold,  we  would  expect 
different  scaling  for  two  reasons.  The  first  is  that  in 
this  case  the  we  instability  lies  in  the  lower  hybrid 
branch  and  is  convective.  The  second  is  the  possiblility 
that  the  beam  radius  as  given  by  (10b)  is  too  short  and 
other  effects  (e.g.  cathode  size,  charge  neutralization 
etc.)  will  control  it.  To  our  surprise  a  recent 
experiment  by  Kawashima22  with  B^l-1.5  kG,  not  only  shows 
scaling  close  to  Equ.  (14),  but  also  the  observed  values 
Ic  are  to  within  less  than  a  factor  of  two  from  the 
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predicted  ones.  The  only  explanation  we  can  provide  is 
the  possibility  that  reflections  from  boundaries  allow  for 
oscillator  behavior  under  conditions  given  by  Equ.  (9). 
This  was  actually  observed  in  a  BPI  experiment  by  the 
Steven's  group. 

Is  it  possible  to  spark  B PD  with  a  convective 
instability  near  ->e ?  The  answer  to  this  question  is  an 
unqualified  yes,  as  long  as  the  system  size  has  many 
growth  lengths  for  sufficient  energy  deposition.  This 
would  be  the  case  for  ionospheric  injection.  In  most  of 

U> 

the  regions  of  interest  to  ionospheric  injection  >1 

’  e 

without  beam  ionization.  Currents  much  smaller  than  in 
the  laboratory  will  be  sufficient  to  trigger  BPD  if  the 
Towsend  condition  is  fulfilled  since  the  system  is  long 
enough  not  to  require  absolute  instability.  To  trigger 
BPD  in  the  laboratory  for  — ~  <1  with  convective 

a 

e 

amplification  would  probably  require  very  strong  growth  in 

n 

the  sense  of  >0(1). 
n 
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FIGURE  CAPTIONS 

Critical  current  for  BPD  ignition  as  a  function, 
of  neutral  gas  pressure  for  three  values  of  the 
magnetic  field  (from  Ref.  7). 

Diagram  of  physical  process  controlling  BPD 
phenomena. 
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FIGURE  2. 


